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1. INTRODUCTION 

Electric motors are generally used in many systems, including household, military, medical and 
industrial [1], [2]. Direct current (DC) motors are considered one of the most important electric motors [3]— 
[5]. They have specifications that make their choice to move many machines within different applications 
[6]—[8]. It is the most appropriate technically with high performance and spatiality, occupying less place due 
to its small size [9]-[11]. DC motors have occupied a wide range of application fields due to their high 
performance. DC motors include conventional DC and brushless DC electric motor (BLDC) [12]-[14]. 

Conventional control systems handle errors resulting from the work of linear systems, any 
adjustment process for their parameters according to the linear system is designed [15]-[17]. An algorithm is 
developed to suit the system that needs a control unit to improve its performance [18]-[20]. In order to 
improve the performance of the work of linear systems and the shortcomings of traditional performance 
[21]-[23]. Many solutions were found, including the advanced optimization [24]-[26]. They simulate the 
hunter and prey and the iteration process as a result of similar cases and experience, including the genetic 
algorithm (GA) and Particle swarm optimization (PSO) as well as the wolves grey wolf optimization (GWO) 
and there are other ways [27]-[29]. 

In the current work, the GWO algorithm was chosen to adjust parameters the traditional 
proportional—integral—derivative (PID) controller to improve the working performance of the DC motor. To 
control the motor speed no it is necessary to specify a target function such as integral time absolute error 
(ITAE). The work of this type of techniques depends on the creation of a random number of the hunting 
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group and the target function determines the hunting behavior and works to estimate the location of the prey 
through iterative processes to obtain the optimal location of the prey. 


2. DC MOTOR AND SIMULATION MODEL 

The electric motor works to produce mechanical energy through the torque arising from the passage 
of an electric current in the coils of the motor [30]—[32]. First, which in turn acts as an influence in the 
magnetic circuit in the motor and thus produces a magnetic driving force. Second, which in turn generates a 
torque to rotate the moving part [33]-[35] which is the rotating part of the electric motor, which produces 
kinetic energy or mechanical [36]—[38]. The electric motor has special specifications that make the selection 
process for any application easy [39]-[41]. Such as rotational speed, number of poles, capacity, current, and 
voltage. according to the manufacturer, and each quantity has a symbol and amount. Table 1 represents the 
specifications of a DC motor that is intended to perform a computer simulation. To represent it 
mathematically first and to find the transformation function which is used in building the simulation model. 
A working algorithm is developed to simulate the operation of the engine first without quantum units. 
Secondly, the engine is run with the traditional control unit and compared with the previous case. Third, the 
simulations are carried out with the optimized advanced technique and the comparison with the previous 
cases. After the simulation, it is possible to work on the appropriate design to control the work of the electric 
motor by setting an algorithm for better performance, higher speed and accuracy in quality to reach a state 
close to ideal in performance. In Figure 1 show the simple model for DC motor. In Figure 2 show the block 
diagram for DC motor at Tf = 0 (no load). In Figure 3 show the block diagram for DC motor with Tf (load). 


Figure 1. Simple model for DC motor 
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Figure 2. Block diagram for DC motor at Tf = 0 (no load) 
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Figure3. Block diagram for DC Motor with Tf(load) 


The electric motor has an electric part and it can be put into a mathematical model with equations 
through which the theoretical specifications can be obtained and the ideal calculations are calculated through 
it. Equations (1) and (2) include both voltage and power, respectively. The electric motor has a mechanical 
second part that can be represented in equations (3) and (4) for speed and torque straight [42]—[44]. 
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ea = İa. Ra — €bV 


ep = Ke. OmV (1) 
Where, K. = back emf constant, wm = motor speed (rad/sec) 


ea = (1.2 * 2) + (0.06 * 500) = 32.4 Volt 
P = Wm. TV (2) 


Where, P = shaft power, T = motor torque 
P = 500 * 0.108 = 54 watt 


Where, ea(t) = armature terminal voltage, i,(t) = armature current, Ra = armature resistance, La = armature 
inductance, ep(t) = back emf. 


kre- QTR 


Om BeK 
Wm = SE 8) 
_ Kyea + (Tp. Ra 
m Boke 


wm = (0.06 * 32.4) + [(1.2 * 0.012)/(0.06 * 0.06)] = 536 rad/s 

Where, Kr = torque constant, T;= static friction torque, jm= rotational inertia, Bm = viscous friction. 
T = Ky.ig — TN -m (4) 
T = (0.06 * 2) — 0.012 = 0.108 N — m 


The mathematical model of the electric motor can also be described by the mathematical equations: 


Electrical equation: e, (t) = Ra. ia(t) + pao + e,(t) (5) 
Lap. Tra. of electrical equation: E,(s) = L.s.I,(s) + Ra + la(s) (6) 
la(s) = pp; [Ea(s) - Ey(S)] (7) 
Mechanical equation: T(t) = jm m O+ Bn: Om (t) (8) 
Lap. Tra. of mechanical equation: T(s) = (jm. Ss + Bm). Wm(S) (9) 
Wm(S) = [1/Gm-S + Bm)]- T(S) (10) 
Electromechanical relationships: e,(t) = Kg. w(t) (11)T(t) = Ky. i, (t) (11) 


Lap. Tra. of electromechanical relationships:e,(s) = Kg. @m(S) (13)T(s) = Ky. ig(s) (12) 


A transfer function for a DC motor as show in equation 13 [43]. 


wm (s) = G(s) = Kr (13) 
Ea(s)  1+G(s).H(s)  Lajm.S2+(Rajm+BmLa).s+(Kr Kg+Ra.Bm) 


16.13 


T.F for DC motor = ——— z 
0.00333s*+0.201s+1 


(14) 
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Table 1. Parameters of DC motor 


Parameters of DC motor Meaning 

Wm 500 rad/sec 

Iq 2A 

Ke 0.06 V.s/red 

Kr 0.06 N.m/A 

Te 0.012 N-m 

Ra 1.2 ohm 

jm 6.2*10-4 N-m-s/rad 
Bn 1*10-4 N-m-s/rad 
La 0.02 H 


Hierarchical algorithm, as shown in Figure 4, is called wolf algorithm and can be represented by 
categories that include the leader, who is from the first level to the top of the hierarchy and is called alpha (a). 
This category is concerned with making decisions such as the spatial environment for living and hunting and 
can be male or female or both. The second category is a body considered in the second level of the hierarchy. 
Hierarchical duty helps the first category to make some decisions and is called beta (B). The third category is 
lower than its predecessors in terms of hierarchy called delta (5), where it is considered a category subject to 
the orders of other groups and has a role to play, which is being a scapegoat. Below the pyramid is the lowest 
rank category, the fourth, which is the category of sheikhs, which performs the functions of guard and scouts. 
It is considered dominant over its predecessor, omega (), but is subject to both alpha and betas. In Figure 4 
show the hierarchy of GWO. Also, the GWO had flow chart that show in Figure 5. 


Alpha 
Beta 
EN » Delta 
Omega oe 


Figure 4. Hierarchy of GWO 
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Figure 5. Flow chart for GWO 


3. SIMULATION MODEL FOR DC MOTOR 
3.1. Simulation model for DC motor with constan speed 

The objective are design of PID controller and tuned it by GWO with transfer function (T.F). The 
T.F of the system in open loop as show in (16). 


16.13 


T.F for DC motor = 
Mie motor = 00033352 + 0.201s + 1 


In Figures 10-12 the simulation model for DC motor. First the simulation model for DC motor 
without controller as show in Figure 10. Second the simulation model for DC motor with PID controller as 
show in Figure 11. Then, simulation model for DC motor with GWO-PID controller as show in Figure 12. 

To improve the performance of the work of systems asks for the search for performance indications, 
which are placed to be indications to improve the performance of systems. To design appropriate controllers 
for any system are tested through one or more performance indications for a high effective system and within 
the required specifications. These indications can be recognized through the following equations for each 
type: i) Integral absolute error (IAE) as show in (15). Also in Figure 6 show block diagram forsimulation 
model of IAE; ii) Integral square error (ISE) as show in (16). Also in Figure 7 show block diagram for 
simulation model of ISE; iii) Integral of the product between the squared time and absolute error (ISTAE) as 
show in (17). Also in Figure 8 show block diagram for simulation model of ISTAE; and iv) ITAE integrates 
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the absolute error multiplied by the time as show in (18). Also in Figure 9 show block diagram for simulation 
model of ITAE. 


IAE = fle(t)|dt (15) 


: 
AY 


Figure 6. Block diagram for simulation model of IAE 


ISE = f e?(t)dt (16) 


Sra iG} 


Figure 7. Block diagram for simulation model of ISE 


ISTAE = f t?]e(t)|dt (17) 
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Figure 8. Block diagram for simulation model of ISTAE 


ITAE = f tle(t)|dt (18) 


Mnene 


Figure 9. Block diagram for simulation model of ITAE 


Figure 10 shows the block diagram for simulation model for DC motor without controller. Figure 11 
shows the block diagram for simulation model for DC motor with PID controller. Figure 12 shows the block 
diagram for simulation model for DC motor with GWO-PID controller. 
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0.003335? +0.201s + 1 


Figure 12. Block diagram for simulation model for DC motor with GWO-PID controller 


3.2. Simulation model for DC motor with variable speed 

In this section there are three parts, in Figures 13-15 simulation model for DC motor of variable 
speed. First, the simulation model for DC motor of variable speed without controller as show in Figure 13. 
Second, the simulation model for DC motor of variable speed with PID controller as show in Figure 14. Last, 
the simulation model for DC motor of variable speed with GWO-PID controller as show in Figure 15. 


GD 16.13 
s? + 0.201s + 0.00333 
DC motor without Controller 


Figures 13. Simulation model for DC motor of variable speed without controller 
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DC motor with PID Controller 
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Figures 14. Simulation model for DC motor of variable speed without controller 
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Figures 15. Simulation model for DC motor of variable speed with GWO-PID controller 


4. 


SIMULATION RESULTS OF DC MOTOR 


4.1. Simulation results of DC motor with constant speed 


This section includes three parts. First by using the simulation model in Figure 10. In this part the 
simulation results as show in Figure 16. Second by using the simulation model in Figure 11. In this part the 
simulation results as show in Figure 17. Then, by using the simulation model in DC motor with GWO-PID 
controller in this part the simulation results as show with comparative among for these simulation pars as 


show in the Table 2 and in Figure 18. 
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Figure 16. Simulation results for DC motor without controller 
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Figure 18. Simulation results for DC motor without controller, with PID and GWO-PID controller 
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Figure 17. Simulation results for DC motor with PID controller 
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Table 2. Simulation results with different algorithm 


Algorithms Over shoot (%) Under shoot (%) Rise time (sec) 
DC motor without Controller 22.840 5.675 0.015 
DC motor with PID Controller 5.851 1.871 0.0021 
DC motor with GWO-pid controller - 0.002 


4.2. Simulation results of DC motor with variable speed 

In this section there are three parts, in Figures 19-21 simulation results for DC motor of variable 
speed. First the simulation response for DC motor of variable speed without controller as show in Figure 19. 
Second the simulation response for DC motor of variable speed with PID controller as show in Figure 20. 
Third the simulation response for DC motor of variable speed with GWO-PID controller as show in 


Figure 21. 
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Simulation response without controller 
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Figure 19. Simulation response for DC motor of variable speed without controller 


Simulation response with PID controller 
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Figure 20. Simulation response for DC motor of variable speed with PID controller 


Simulation response with GWO-PID 
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Figure 21. Simulation response for DC motor of variable speed with GWO-PID controller 


In Figure 22 show the ITAE with each iteration. In Figure 23 show the step response. In Figure 24 
show the wondow of T.F code for DC motor. In Figure 25 show the wondow of code for GWO-PID. 
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Figure 23. Step response 
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Figure 24. Window of T.F code for DC motor 
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Figure 25. Window of code for GWO-PID 


5. CONCLUSION 

The current study was conducted by simulating the operation of a DC motor with working 
conditions for three cases that were selected and suggested to verify the preference for using the optimal 
advance and compare it with the two cases of no control systems. Others with the presence of control using 
the traditional system, a traditional control unit. The results demonstrated the superiority of the optimum 
advanced over the traditional progress unit in terms of response speed, time to reach a stable state, accuracy 
and quality. To adjust the parameters of the traditional controllers using the optimum advanced, an 
appropriate mechanism and technology from the advanced optimization techniques were chosen, as the Gray 
wolf technology algorithm was chosen as an optimization technique and ITAE to adjust the parameters of the 
traditional PID controller. The response speed, high accuracy, stability time and overtaking are higher and 
lower compared to different cases of working conditions that simulate real time. 
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